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ABSTRACT 



Aims. We gathered about 100 high-resolution spectra of three typical HgMn (mercury-manganese) stars, HD 11753, HD 53244, and 
HD 221507, to search for slowly pulsating B-like pulsations and surface inhomogeneous distribution of various chemical elements. 
Methods. Classical frequency analysis methods were used to detect line profile variability and to determine the variation period. 
Doppler imaging reconstruction was performed to obtain abundance maps of chemical elements on the stellar surface. 
Results. For HD 1 1753, which is the star with the most pronounced variability, distinct spectral line profile changes were detected for 
Ti, Sr, Y, Zr, and Hg, whereas for HD 53244 and HD 221507 the most variable line profiles belong to the elements Hg and Y, respec- 
tively. We derived rotation periods for all three stars from the variations of radial velocities and equivalent widths of spectral lines 
belonging to inhomogeneously distributed elements: P ro , (HD 1 1753)=9.54 d, P ro , (HD 53244)=6.16 d, and P ro( (HD221507)=1.93 d. 
For HD 1 1753 the Doppler imaging technique was applied to derive the distribution of the most variable elements Ti, Sr, and Y using 
two datasets separated by ~65 days. Results of Doppler imaging reconstruction revealed noticeable changes in the surface distribu- 
tions of Ti a, Sr n, and Y n between the datasets, indicating the hitherto not well understood physical processes in stars with radiative 
envelopes that cause a rather fast dynamical chemical spot evolution. 
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1. Introduction 

The mercury-manganese (HgMn) stars constitute a well-defined 
sub-class of chemically peculiar (CP) stars of the B7-B9 spectral 
types with T e g between 10000 and 15 000 K. These stars exhibit 
marked abundance anomalies of several elements: e.g., over- 
abundances of Hg, Mn, Ga, Y, Cu, Be, P, Bi, Sr, Zr, and deficien- 
cies of He, Al, Zn, Ni, Co. More than two thirds of them belong 
to spectroscopic binaries (Hubrig & Mathys 1995). They are 
slow rotators ((v sin/) = 29 kms -1 , Abt | 19721 . There is no ev- 
idence that they would have strong large-scaled organized mag- 
netic fields. Their elemental overabundances/underabundances 
are believed to be due to radiatively-driven diffusion and gravi- 
tational settling. 

In the H-R diagram many HgMn stars are located in the in- 
stability strip of the so-called slowly pulsating B (SPB) stars (see 
De Cat 2003 for a review on the latter), and the sophisticated 
models predict that pulsations should also be driven in HgMn 
stars (Turcotte & Richard 2005). Searches for variability in this 
group of stars have been made mostly photometrically in the 
past, but without any success. Very recently, Alecian et al. (2009) 
discovered low amplitude (less than 1 .6 mmag) periodic varia- 
tions (4.3 and 2.53 days respectively, with harmonics) in two 
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candidate HgMn stars by means of the high quality light curves 
provided by the CoRoT satellite. These variations are compati- 
ble with theoretically predicted pulsation periods. However, as 
stated by the authors, only spectroscopic datasets could help to 
conclusively establish or withdraw this pulsation interpretation. 
Currently, there is thus no observational proof of pulsation in 
HgMn stars. 



The aspect of inhomogeneous distribution of some chemi- 
cal elements over the surface of HgMn stars was for the first 
time discussed by Hubrig & Mathys (1995). From the survey 
of HgMn stars in close spectroscopic binaries (SBs) it was sug- 
gested that some chemical elements might be inhomogeneously 
distributed on the surface, with in particular a preferential con- 
centration of Hg along the equator. In close double-lined systems 
(SB2s), where the orbital plane has a small inclination to the line 
of sight, a rather large overabundance of Hg was found. By con- 
trast, in stars with orbits almost perpendicular to the line of sight, 
mercury is not observed at all. The first indication of variabil- 
ity of the Hgll 3984 A and YII 3983 A lines wa s repo rted for 
the HgMn SB2 system ARAur by Takeda et al. dl9791 l. Later, 
Wahlgren et al. ( 1200 11 1 and Adelman et al. (2002]) showed that 
the Hg II 3984 A line of the primary component of a And varies 
with a 2.8-d period. The spectral line variations were attributed 
to the surface inhomogeneous mercury distribution along the 
stellar equator, together with the stellar rotation period. 
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Table 1. Atmospheric parameters (effective temperature, loga- 
rithm of the gravity, microturbulent velocity, and projected rota- 
tional velocity) taken from Dolk et al. (20031 for HD 1 1753 and 
HD 221507, and from Woolf & Lambert (1999 ) for HD 53244. 
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Position of the studied stars in a (log T e g Jog g)-diagram. 
left to right, the error boxes represent HD 53244, 
HD 221507, and HD 1 1753. The full and dashed lines represent 
the boundary of the theoretical SPB instability strip for a metal- 
licity Z = 0.02 and Z = 0.01, respectively (taken from Miglio et 
al. 120071) . 
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Fig. 2. Spots on the surface of HD 11753: time-series spectra 
phased on the stellar rotation period of 9.54 days. 
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Fig. 3. Left: Inhomogeneous distribution of Hgn on the sur- 
face of HD 53244 apparent in the time-series spectra around 
Hgn /13984 phased on the rotation period of 6.16 days. Right: 
Inhomogeneous distribution of Y n on the surface of HD 22 1 507 
apparent in time-series spectra around the spectral line of Y n 
/I4310 phased on the rotation period of 1.93 days. 



Recently, Kochukhov et al. ((2005 ) found clear signatures of 
surface mercury spots in two rapidly rotating HgMn stars by 
analysing the Hgll 3984 A line profiles. Variability of spectral 
lines associated to larger number of chemical elements were dis- 
covered for the first time by Hubrig et al. (2006a) for the primary 
component of the eclipsing binary AR Aur. The strongest vari- 
ations were found for the chemical elements Pt, Hg, Sr, Y, Zr, 
He, and Nd. The first Doppler maps for the elements Mn, Sr, Y, 
and Hg were recently presented by Savanov et al. (2009). The 
study of Hubrig et al. (2008 ) suggests that spectral variability 
of various chemical elements is indeed observed in most HgMn 
stars. 

We present the first observational study based on a substan- 
tial number of spectra, more than one hundred, obtained with the 
CORALIE echelle spectrograph attached to the 1.2m Leonard 
Euler telescope in La Silla in Chile. The selected targets, the 
single-lined (SB1) spectroscopic binaries HD 11753 (0Phe, V 
= 5.1 mag, B9p) and HD 53244 (yCMa, V = 4.1 mag, B8II), 
and the star HD 221507 (J3 Scl, V = 4.4 mag, B9.5IVmnpe) 
were chosen as the brightest known southern HgMn stars visi- 
ble during the periods of observation. The goal of the presented 
spectroscopic study was twofold: to search for stellar pulsations 
and/or surface inhomogeneous distribution of various chemical 
elements. Our observations and data reduction are presented in 
the appendix A. 



2. Line profile variability 

Atmospheric parameters and abundances of several elements 
of selected HgMn stars were studied in the past by Smith & 
Dwore tsky JTf93l ), Woolf & Lambert ([1999] ) and Dolk et al. 
(2003). The stellar parameters for our three studied stars are 
given in Table [TJ Their position in a (log T e g Jog g)-diagram is 
presented in Fig.Q] The star HD 221507 is located in the SPB 
instability strip, whereas HD 11753 and HD 53244 are located 
outside this strip. According to theoretical models (e.g. Miglio et 
al. 120071 ) HD 221507 is expected to exhibit SPB-like pulsations. 
To search for variability we studied the behaviour of the radial 
velocities and equivalent widths of the spectral lines of various 
elements. We used the Lomb-Scargle method (Scargle 1982 ) and 
the phase dispersion minimization method (Stellingwerf |1978l ). 
Additionally, we performed a two-dimensional frequency anal- 
ysis, which is available in FAMIAS (Frequency Analysis and 
Mode Identification for Asteroseismology) (Zima 2008 ). 

Among the studied stars, the spectra of HD 11753 exhibit 
the most prominent variability. This was most clearly detectable 
for the elements Ti, Sr, and Y. The Fe lines, by contrast, show 
no profile variations. For that reason, we determined the stellar 
radial velocity (RV) by measuring 16 Fe lines free of blends. 
The star HD 11753 is known as a single-lined binary star, but 
its orbital parameters are not well known. According to our ob- 
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SET2, Till 441 7,71 



SET2, Yll 4900.1200 




Fig. 4. Ti abundance map of HD 11753 obtained from Tin line 
4417.7136A for setl and set2. The colour indicates the abun- 
dance with respect to the total number density of atoms and ions. 



servations the orbital period would be long. In fact, all our RV 
measurements, which extend over 10 months, range from 14.1 to 
14.5 km s _1 . In particular for the two first datasets (76 and 28 ob- 
servations, respectively) the mean RVs are 14.24 + 0.01 kms -1 
(RMS (root mean square)=0.05kms _1 ) and 14.36 + 0.01 kms 1 
(RMS=0.05 kms 1 ). Considering the small dispersion and the 
absence of any trend with the rotational velocity found from Y 
lines, we adopted these values as the stellar RV for the obser- 
vations of the corresponding datasets. For HD 53244 our RVs 
present an RMS of 0.2 km s without significant differences be- 
tween different runs. Therefore we did not apply an RV correc- 
tion other than the subtraction of the mean RV. 

For HD 11753, the radial velocities and equivalent widths 
of the Ti, Sr, and Y lines were found to vary with the period 
P=9.54d. In the spectra of HD 53244 the variations are appar- 
ent for the elements Hg and Mn with the P=6.16 d, while in the 
spectra of HD 221507 variations associated to Hg, Mn, and Y 
are detected, indicating the period P= 1 .93 d. 

For all three observed stars the behaviour of the line pro- 
files is different for different elements. Moreover, spectral lines 
associated to certain elements, like silicon, are constant. These 
characteristics do not support the SPB-like pulsation interpreta- 
tion (De Cat 2001). We thus conclude that the determined pe- 
riods correspond to rotation periods of studied HgMn stars and 
the detected line profile variability is caused by inhomogeneous 
distribution of chemical elements on their stellar surface. Apart 
from two other HgMn stars, a And and AR Aur, no other rotation 
period determinations for HgMn stars were published before this 
study. In Figs.|2]and[3]we display time-series spectra phased on 
the stellar rotation periods around selected lines of various ele- 
ments in all three stars. The images were created by averaging 
all the spectra within the phase interval <p ± 0.04 for each phase 
</>■ 

3. Surface chemical inhomogeneities by Doppler 
imaging reconstruction 

For the star HD 1 1753, which exhibits the most distinct spectral 
line profile variations for several elements, we used the Doppler 
imaging technique to reconstruct the surface distributions of Ti, 
Sr, and Y. The elements Ti and Y have numerous transitions in 



Fig. 5. As in Fig.g] but now using the Yn line 4900. 1200 A. 



SET1, Srll 1215.7190 




SET2, Srll 4215.5190 




Fig. 6. As in Fig.|4] but now using the Srn line 4215.5190 A. 



the observed optical spectral region allowing us to select un- 
blended spectral lines which show strong variability. 

The two sets of observations of HD 1 1753 obtained in 2000 
Sep 28 - Oct 1 1 (setl) and in 2000 Dec 02 - Dec 15 (set2) con- 
sist of 76 and 28 observations, respectively, evenly spread over 
the stellar rotation cycle. Using two separate surface abundance 
maps based on observations of two data sets that are on average 
65 days apart allowed us to obtain important information on the 
temporal evolution of elemental surface inhomogeneities. The 
Doppler imaging technique takes advantage of the partial reso- 
lution of the stellar surface provided by the rotational Doppler 
effect, inverting a line profile time series into a 2-D map of the 
stellar surface. The inversion of time series of spectroscopic ob- 
servations is based on regularised image-reconstruction proce- 
dures implemented in the Doppler imaging code INVERS7PD 
written by N. Piskunov (see, e.g., Piskunov et al. 1990) and mod- 
ified by T. Hackman d2DTHV 

The observations were compared to a grid of local line 
profiles calculated with the SPECTRUM spectral synthesis 
code (Gray & Corbally 1994) and Kurucz model atmospheres 
(Kurucz 1993). The local line profiles were calculated for 10 
limb angles, and the stellar parameters were fixed to the val- 
ues given in Table Q] For inclination, vsin i, and microturbu- 
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lence inversion using several different values was carried out 
to determine the value that best fitted the observations. As best 
values for Ti and Y surface distribution reconstruction we ob- 
tained v sin /=13.5 km s _1 , z'=53°, and £=0.5 km s _1 . On the other 
hand, the best fit for the reconstruction of the Sr distribution was 
achieved for v sin i= 12.3 km s _1 . The discrepancies in the v sin i 
values are very likely related to the Sr vertical abundance stratifi- 
cation, which is frequently observed in chemically peculiar stars 
(e.g. Kochukhov et al. 12006). For each studied element a grid of 
abundances spanning from -4.0 to -9.5 was calculated, using a 
step of 0.5 in abundance. These abundances are in the scale used 
by SPECTRUM and are thus expressed with respect to the total 
number density of atoms and ions, and not with respect to hy- 
drogen, with the log of the abundance of hydrogen set equal to 
12.0. 

The Ti n line chosen for the inversions is 4417.7136 A, which 
has an excitation potential of 1 .237 eV. Chemical maps were sep- 
arately recovered for both datasets (setl and set2). As can be 
seen in Fig. |4j the Ti abundance does not exhibit a distinct ring 
structure around the stellar equator as was found for the Hg dis- 
tribution in the HgMn star aAnd (Adelman et al. ((2002). The 
maps in general show a surface abundance that is higher than 
the solar abundance of Ti, -7.02, and have an average abundance 
of -6.47. Two main structures are well noticeable on the surface: 
a high abundance spot at high latitudes (40° - 75°) spanning 
the phases 0.75-1.00, with an extension towards and beyond the 
phase 0.5 forming a half ring, and a lower abundance spot in the 
polar regions at the phases 0.1-0.5. Also, some discrete spots 
of lower and higher abundance appear in the equatorial region. 
The Ti lower abundance spot at phases 0.1-0.5 is less promi- 
nent in set2 than in setl. Also the spot configuration of the high 
abundance high latitude half ring of Ti at the latitudes 40° - 75° 
changes between the two sets. 

The maps obtained from the Yn line 4900. 1200 A with an 
excitation potential of 1.033 eV is presented in Fig. [5] The maps 
show a high abundance region at phases 0.5-1.0 extending from 
the latitude 45° to the pole, with an extension to the equator 
around phase 0.8. The Y abundance distribution shows a high 
latitude lower abundance spot around phases 0.2-0.4, similarly 
to the Ti abundance maps. Some lower and higher abundance 
spots are also seen at the equatorial region. The average abun- 
dance of the Y maps is -7.01, which is significantly higher than 
the solar abundance of -9.80. We note that all the features re- 
vealed in the maps show abundances that are higher than the 
solar abundance of Y. Similar to the Ti maps, we observed in 
the Y maps that the lower abundance high latitude feature at 
phases 0.2-0.4 becomes more extended and less prominent in 
set2. Additionally, the abundance of the high abundance spot at 
phases 0.6-1.0 increases. 

The Sr maps based on the modelling of the Srn4215.5190A 
line reveal similarities in the distribution of overabundant and 
underabundant features, though the similarity is more pro- 
nounced in the Sr and Y maps. In Fig.|6]the equatorial and polar 
belt-like Sr structures resemble rather well those of Y. The time 
evolution of the chemical spots is also seen in the Sr maps, with 
the lower abundance high latitude feature at phases 0.1-0.5 be- 
coming more extended in set2, and the high abundance feature of 
the phases 0.8-1.0 moving towards the equator and phases 0.5. 
The average abundance of the Sr maps is -6.89, which exceeds 
the solar abundance of -9.07. 

One has to keep in mind though that set2 has a gap at the 
phases 0.187-0.391. This phase gap can affect the achieved reso- 
lution on the surface at these phases, and can especially hamper 



the exact determination of the latitudes. Still, the changes be- 
tween the two sets are seen from the maps at other phases, and 
also from the line profiles themselves as shown in AppendixICl 



4. Discussion 

All Ti, Sr, and Y abundance maps reveal a structure reminiscent 
of broken rings of low and high abundance. This elemental dis- 
tribution is to some extent similar to the maps previously recon- 
structed for another HgMn star, ARAur (Savanov et al. [2009), 
where the elements Mn, Y, Sr, and Hg show abundance con- 
centration in equatorial and polar features. Typically, inhomoge- 
neous chemical abundance distributions are observed only on the 
surface of magnetic chemically peculiar stars with large-scale 
organised magnetic fields. In these stars, the abundance distri- 
bution of certain elements is non-uniform and non-symmetric 
with respect to the rotation axis. A magnetic field of the or- 
der of a few hundred Gauss was detected in hydrogen lines of 
four HgMn stars by Hubrig et al. (2006b) using low-resolution 
(R = 2000) circular polarisation spectra obtained with FORS 1 
at the VLT. This small sample of HgMn stars also included the 
spectrum variable HgMn star a And with a magnetic field of 
the order of a few hundred Gauss. On the other hand, high- 
resolution spectropolarimetric spectra of some HgMn stars, in- 
cluding a And, were used in studies of Shorlin et al. (2002) and 
Wade et al. (2006), where no detection was achieved using all 
metal lines together in the least-squares deconvolution multi- 
line profile. Although strong large-scale magnetic fields have not 
generally been found in HgMn stars, it has never been ruled out 
that these stars might have tangled magnetic fields of the order of 
a few thousand Gauss with no net longitudinal component (e.g., 
Mathy s & Hubrig [1995] Hubrig et al. [19991 Hubrig & Castelli 
1200 II ). It is of interest that magnetohydrodynamical simulations 
by Arlt et al. (2003 ), which combine a poloidal magnetic field 
and differential rotation can produce a magnetic field topology 
that is similar to the broken elemental ring structures seen in 
HD 1 1753 and ARAur. These simulations and their implication 
have been recently discussed by Hubrig et al. ( 120081 ). 

The abundance maps of HD 11753 presented in this work 
exhibit clear differences between the surface abundance distri- 
bution of Ti, Sr, and Y. We also detected clear differences in the 
spot configurations obtained from the same lines but for differ- 
ent data sets, which indicates a rather fast dynamical evolution 
of the abundance distribution with time. Kochukhov et al. (2007 ) 
discovered mercury clouds in the atmosphere of a HgMn star 
a And that showed secular changes with a time period of 2-4 
years. In our analysis, using two datasets separated by ~65 days, 
we reveal that the changes in the chemical spot configuration of 
HD 11753 appear much faster and can already be detected at a 
time scale of months. 

The results reported in this paper open up new perspectives 
for our knowledge and understanding of HgMn stars. Different 
dynamical processes take place in stellar radiation zones. An 
interaction between the differential rotation, the magnetic field, 
and the meridional circulation could possibly play a role in the 
generation of dynamical evolution of chemical spots. From the 
comparison of maps we find that it is possible that the Y and 
Sr distributions show indications of an increasing rotation rate 
towards the rotation pole, so-called differential rotation of anti- 
solar type. On the other hand, further analyses of the elemental 
surface distribution in a larger sample of HgMn stars should be 
carried out before the implication of these new results can be 
discussed in more detail. 
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Appendix A: Observations and data reduction 

Our spectroscopic data were obtained with the CORALIE 
echelle spectrograph attached to the 1.2m Leonard Euler tele- 
scope in La Silla in Chile. In total, we gathered 113, 103, and 
105 useful spectra for HD 11753, HD 53244, and HD 221507, 
respectively. The observing logbook is given in Table I A. 1 1 and 
the obtained phases for HD 1 1753 in Fig. I A. 1 1 The wavelength 
domain of the CORALIE fibre-fed spectrograph ranges from 
3875 to 6820 A recorded on 68 orders. The CCD camera is a 2k 
x 2k CCD with pixels of 15 fim. CORALIE reaches a spectral 
resolution of 50000 with a 3 pixel sampling. 

To search for variability due to pulsation in line profiles 
we used rather short integration times to obtain a S/N ratio of 
about 100 near the wavelength 4130 A, which is the position of 
a prominently present Si II doublet. This spectral region is well- 
suited for studying line profile variability of SPB stars (e.g., De 
Ridder etal. l2U02l . 

We used the online reduction package available for the 
CORALIE spectrograph based on the method of Baranne et al. 
( 119961 1. The process involves the usual steps of de-biasing, flat- 
fielding, background subtraction and wavelength calibration by 



Table A.l. Logbook of our spectroscopic observations of the 
three HgMn stars. N denotes the number of observations and the 
HJD is given with respect to HJD = 2400000. 
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Fig. A.l. Phases of the HD 11753 observations for both sets (setl 
on the left). 



means of measurements of a ThAr calibration lamp. We did a 
more precise correction for the pixel-to-pixel sensitivity varia- 
tions by using all available flat fields obtained during the night 
instead of using only one flat field, as is done by the on-line re- 
duction procedure. After a correction to the heliocentric frame, 
the spectra were normalised using a cubic spline fit. 

Appendix B: Details on Doppler imaging 

The average S/N ratio of the HD 11753 observations for setl 
and set2, measured at 4400 A are 132 and 143, respectively. 
This means that in Doppler imaging a smaller deviation than 
this should not be achieved between the model and the obser- 
vations. In the Tikhonov regularisation used by INVERS7PD 
the regularisation has to be chosen in a way that the noise 
in the observations is not mapped. For HD 11753 the devia- 
tion between the model and observations for the setl is be- 
tween 0.750% (Yn4900.120A) and 0.893% (Srn4215.5190 A), 
and for the set2 between 0.628% (Yn 4900. 1200 A) and 0.803% 
(Srn 4215. 5190 A). These values are comparable to the average 
noise level in the observations. One has to also keep in mind that 
the noise level in the blue part is higher than in the red part of 
the spectrum. Furthermore, the comparison of the model and ob- 
servations shows that the line profiles are well fitted without the 
model following the noise. An example of the line profile fits is 
shown in Fig. |BT]for the Y n 4900. 1200A line of set2. 

Appendix C: Line variability 

The changes between the two sets are also seen when looking at 
the line profiles at the phases where the largest differences in the 




Fig. B.l. Set2 Yn 4900. 1200 A spectra of HD 1 1753 used in the Doppler imaging. Calculated and observed spectral lines are shown 
by lines and crosses, respectively. The numbers on the right side of the plots give the phase for each spectrum. 
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Fig. C.l. Mean YII line profiles from close-by phases for setl (phases 0.176-0.177) and set2 (phases 0.183-0.187). The two mean 
spectra are plotted for three YII lines, which have similar excitation energies: 4900.1200 A, 5087.4160 A and YII 5205.7240 A. In 
the plots setl is given by the solid line and set2 by the dashed line. 



maps occur. Figure IC.ll shows averaged spectra of very close- 
by phases from setl (mean of two spectra, phases: 0.176-0.177) 
and set2 (mean of three spectra, phases: 0.183-0.187). The mean 
profiles for the two sets are shown for three Y II lines of similar 
strength and excitation energy (4900.1200 A, 5087.4160 A and 
5205.7240 A). As can be seen all three lines show changes be- 
tween the two sets that are similar and in-line with the changes 
seen in the Doppler images, i.e., the Y II abundance is higher for 
set2 than setl at these phases. 



